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ABSTRACT 

We measure the evolution of the galaxy stellar mass function from z = 1.3 to z = 0.5 using the first 53 608 redshifts of the ongoing 
VIMOS Public Extragalactic Survey (VIPERS). Thanks to its large volume and depth, VIPERS provides a detailed picture of the 
galaxy distribution at z — 0.8, when the Universe was =i 7Gyr old. We carefully estimate the uncertainties and systematic effects 
associated with the SED fitting procedure used to derive galaxy stellar masses. We estimate the galaxy stellar mass function at several 
epochs between z = 0.5 and 1.3, discussing in detail the amount of cosmic variance affecting our estimate. We find that Poisson 
noise and cosmic variance of the galaxy mass function in the VIPERS survey are comparable with the statistical uncertainties of large 
surveys in the local universe. VIPERS data allow us to determine with unprecedented accuracy the high-mass tail of the galaxy stellar 
mass function, which includes a significant number of galaxies that are usually too rare to detect with any of the past spectroscopic 
surveys. At the epochs sampled by VIPERS, massive galaxies had already assembled most of their stellar mass. We compare our 
results with both previous observations and theoretical models. We apply a photometric classification in the (U - V) rest-frame 
colour to compute the mass function of blue and red galaxies, and we find evidence for the evolution of their contribution to the total 
number density budget: the transition mass above which red galaxies dominate is found to be about 10 104 Me at z — 0.55 and evolves 
proportional to (1 +z) 3 . We are able to trace separately the evolution of the number density of blue and red galaxies with masses above 
10 11 4 Me, in a mass range barely studied in previous work due to the small volumes surveyed. We find that for such large masses, red 
galaxies show a milder evolution with redshift, when compared to objects at lower masses. At the same time, we detect a population 
of similarly massive blue galaxies, which are no longer detectable below z = 0.7. These results show the improved statistical power 
of VIPERS data, and give initial promising indications of mass-dependent quenching of galaxies at z — 1. 

Key words. Galaxies: fundamental parameters, mass function, evolution, statistics - Cosmology: observations 

1. Introduction large amounts of data, thus facilitating a number of important 

statistical studies. The galaxy stellar mass function (GSMF), 
The last decade has seen many significant advances in the study defined afj the co . movina number density of galaxies within a 

k£j of galaxy evolution prompted by large astronomical surveys. In stellar mass bin (M M + dMl h one such fundamen tal statis- 

_ c3 particular, such surveys sample large cosmic volumes collecting ^ a ii owing us t0 trace the history of baryonic mass assembly. 

Measurements of the GSMF help constraining the cosmic star 
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vations obtained with MegaPrime/MegaCam, a joint project of CFHT In the nea rby universe, the GSMF has been measured to 

and CEA/DAPNIA, at the Canada-France-Hawaii Telescope (CFHT), M h acc b exploiting the Two Micron All Sky Survey 

which is operated by the National Research Council (NRC) of Canada (2MASS) thg 2dF Gal Redshift s (2dFGRS, Cole et al. 
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APIX and the Canadian Astronomy Data Centre as part of the Canada- (1976) function, with an upturn at M - 10 W M B (Baldry et al. 

France-Hawaii Telescope Legacy Survey, a collaborative project of 2008; Li & White 2009; Baldry et al. 2012). Such bimodal- 
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2005), is thought to be related to the existence of two distinct 
galaxy types: a population of star-forming galaxies, with blue 
colours and disc-dominated or irregular morphology, and a class 
of red early-type galaxies that, on the contrary, have their star 
formation substantially shut off (Kauffmann et al. 2003a; Franx 
et al. 2008; Bell et al. 2007). 

At higher redshift, such studies are more challenging be- 
cause of the faintness of the objects. However, early seminal 
work took advantage of the Hubble Space Telescope to con- 
struct samples of a few hundreds galaxies up to z — 3, find- 
ing evidence for an increase of the average stellar mass density 
with cosmic time (Rudnick et al. 2003; Dickinson et al. 2003; 
Fontana et al. 2003). Later, deeper surveys were able to show the 
lack of evolution in the high-mass end of the GSMF (GOODS- 
MUSIC catalogue, Fontana et al. 2006) opposed to an increase 
of galaxy density at lower masses (VVDS survey, Pozzetti et al. 
2007), a result consolidated up to z — 4 by means of near- and 
mid-infrared data which facilitate better estimates of the stel- 
lar masses (Perez-Gonzalez et al. 2008; Kajisawa et al. 2009). 
Although some disagreements exist, such findings indicate that 
massive galaxies were assembled earlier than those with lower 
stellar mass, suggesting a 'downsizing in stellar mass' has taken 
place (Fontanot et al. 2009). 

Besides these results, first attempts to study the GSMF by 
dividing blue/active from red/quiescent objects provided inter- 
esting results, despite the relatively limited statistics, and re- 
vealed that within the GSMF the number of blue galaxies at in- 
termediate masses decreases as a function of cosmic time, while 
the fraction of red galaxies increases (Bundy et al. 2006; Borch 
et al. 2006). This early work was extended using larger galaxy 
samples (as in COSMOS and zCOSMOS, Drory et al. 2009; Li- 
bert et al. 2010; Pozzetti et al. 2010), or very deep observations 
(GOODS-NICMOS survey, Mortlock et al. 201 1 ), producing ro- 
bust results for the evolution in number density of both these 
galaxy populations. They also showed that a double Schechter 
function is a good fit to the GSMF data out to to z — 1 (Pozzetti 
et al. 2010; Peng et al. 2010). 

A fundamental picture emerging from these studies is the 
transformation of star-forming galaxies into "red and dead" ob- 
jects through some quenching mechanism that halts the produc- 
tion of new stars. To discriminate between the various quench- 
ing mechanisms proposed in the literature (see Gabor et al. 2010, 
and reference therein), it is crucial to obtain precise and accurate 
measurements to constrain theoretical models (Lu et al. 2012; 
Mutch et al. 2013; Wang et al. 2013). Unfortunately, such com- 
parisons are hard, as on one side the modelling of galaxy evo- 
lution, based on A^-body dark matter simulations (e.g. De Lucia 
& Blaizot 2007; Bower et al. 2006; Guo et al. 2011, 2013), re- 
quires a high level of complexity to parametrise all the physical 
processes (star formation, supernova ejecta, etc.). On the obser- 
vational side, instead, it is hard to derive physical galaxy prop- 
erties with sufficient precision to constrain models, especially 
for rare populations like the most massive galaxies (Marchesini 
et al. 2009, 2010). 

The latest galaxy surveys are helping with accurate measure- 
ments of the GSMF, especially in the high-mass end, and could 
solve discrepancies between data and models (BOSS, Maraston 
et al. 2012). State-of-the-art analyses provide new evidence con- 
firming the dependence on stellar mass of the physical processes 
that extinguish star formation: from z - 1 to z — 3, quenching 
is more efficient for M > 10 10 8 Mq (Ilbert et al. 2013, using 
UltraVISTA data), while at z < 1 the number density of quies- 
cent galaxies evolves significantly at lower masses (Ilbert et al. 
2010; Moustakas et al. 2013 using PRIMUS data). The role en- 



vironment plays in quenching galaxies is still debated (Cucciati 
et al. 2010; Iovino et al. 2010; Bolzonella et al. 2010; Peng et al. 
2010;Vulcani et al. 2013). 

Within this context, the VIMOS Public Extragalactic Red- 
shift Survey (VIPERS) provides a novel opportunity. As we shall 
describe here, this survey provides a combination of wide angle 
coverage, depth and sampling that proves to be ideal to measure 
the GMSF at z ~ 1 with unprecedented precision. The large vol- 
ume allows to effectively probe the massive end of the GSMF 
at these redshifts; compared to photometric redshift surveys, the 
accurate spectroscopic redshift measurements are also crucial to 
for an accurate measurement at the high-mass end, where a few 
interlopers can dramatically change the shape of the GSMF. 

In this paper we present the first measurements of the GSMF 
from the catalogue of the first ~ 55 000 objects; in this first 
analysis we shall concentrate on discussing the evolution of the 
GSMF from z = 1.3 down to z = 0.5, i.e. within the range cov- 
ered by the VIPERS data, for the whole galaxy sample and sep- 
arately for the blue and red populations. We shall also discuss in 
detail the sources of error and potential systematic effects, which 
become dominant at the level of precision on the GSMF allowed 
by the VIPERS data. These are particularly crucial when com- 
parison with independent data is performed. 

In Sect. 2 we present the the VIPERS galaxy catalogue, 
which has been used in this work, and describe how stellar 
masses have been estimated through the SED fitting technique. 
The global mass function is presented in Sect. 3, along with a 
discussion on the sample completeness and the main sources 
of uncertainties. We compare those results with both previ- 
ous surveys and models in Sect. 4. In Sect. 5, after applying 
a colour classification, we study the mass function (and the re- 
lated number density) of red and blue galaxies. Our results are 
summarised in Sect. 6. Unless specified otherwise, our cos- 
mological framework assumes Q„, = 0.25, Qa = 0.75, and 
h 70 = // /(70kms _I Mpc~'). All the magnitudes are in the AB 
system (Oke 1974). 



2. Data 

VIPERS 1 is an ongoing redshift survey aiming at observing ap- 
proximately 100000 galaxies and AGNs at intermediate red- 
shifts ((z) ~ 0.8) in the magnitude range of 17.5 < ; < 
22.5. At the completion of the survey, expected in 2014, ap- 
proximately 24deg 2 will have been covered within two fields 
of the Canada-France-Hawaii Telescope Legacy Survey Wide 
(CFHTLS-Wide) 2 , namely Wl and W4. The sky region covered 
at present is ~ 7.5 deg 2 in each of Wl and W4, with an effec- 
tive area of 5.34 deg 2 and 4.97 deg 2 respectively after account- 
ing for the photometric and spectroscopic masks, or ~ 10 deg 2 
as a whole. Once completed, VIPERS will be the largest spec- 
troscopic survey at such redshifts in terms of volume explored 
(1.5 x 10 8 Mpc /iyo). All details on the survey design and con- 
struction can be found in Guzzo et al. (2013). 

The main science drivers of VIPERS are the accurate mea- 
surement of galaxy clustering, bias parameter, and the growth 
rate of structures, along with the study of the statistical proper- 
ties of galaxies and their evolution when the Universe was about 
half its current age. These topics are the subject of the paral- 
lel accompanying papers of this series (Guzzo et al. 2013; de la 
Torre et al. 2013; Marulli et al. 2013; Fritz et al. 2013; Malek 



1 http://vipers.inaf.it 

2 http : //www . cfht . hawaii . edu/Science/CFHLS/ 
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et al. 2013; Bel et al. 2013). A previous smaller VIPERS sam- 
ple has already been used in Granett et al. (2012) and Marchetti 
et al. (2013). 

The spectroscopic survey is complemented by a set of pho- 
tometric ancillary data, obtained from public surveys and dedi- 
cated observations, allowing us to estimate with high precision 
several galaxy properties, in particular galaxy stellar masses and 
rest-frame magnitudes. 

2.1. Photometry 

The VIPERS spectroscopic sample has been selected from the 
Wl and W4 fields of the CFHTLS-Wide. Therefore, for each 
galaxy we have a photometric dataset consisting of u*, g' , r', 
/', and z' magnitudes (SExtractor's MAG_AUTO derived in 
double image mode in order to maintain the same aperture in all 
bands, Bertin & Arnouts 1996), as measured by the Terapix team 
for the T0005 data release (Mellier et al. 2008). The Terapix 
photometric masks, which discard areas around bright stars or 
with problematic observations, have been revisited by our team 
to recover regions where the photometric quality is deemed suf- 
ficient for our analysis (Guzzo et al. 2013). 

We took advantage of the full wavelength range of the 
VIPERS photometric dataset, as this significantly improves the 
results of our SED fitting; in particular, near-infrared (NIR) 
fluxes are critical to constrain physical parameters and break de- 
generacies between the mean age of the stellar population and 
dust attenuation, and they allow to compute robust estimate of 
stellar masses (e.g. Lee et al. 2009). 

In order to exploit the full potential of VIPERS in analysing 
the galaxy properties as a function of time and environment, we 
have undertaken a follow-up in the /if -band in the two VIPERS 
fields with the WIRCAM instrument at CFHT and in the far- 
and near-UV (FUV and NUV) channel with the GALEX satel- 
lite (Arnouts et al., in prep.). The /if -band observations were 
collected in between 2010 and 2012 with several PI and discre- 
tionary time programs. The /if-band depth has been optimised to 
match the brightness of the spectroscopic sources: at the mag- 
nitude limit (K wircam — 22.0 at 5<x), 96% of the spectroscopic 
sample in W4 is observed in K wircam, while in Wl this percent- 
age is approximately 80%. 

In addition to WIRCAM data, we match our CFHTLS op- 
tical catalogue with the recent UKIDSS data releases 3 using a 
matching radius of 0.8". The Wl field overlaps with UDS and 
DXS, whereas the W4 field is fully covered by the shallower 
LAS and partially covered by DXS. We use, where available, 
Petrosian magnitudes in the Y, J, H and K bands converted in 
the AB system. Considering also ZifuKiDSS, the percentage of our 
spectroscopic sample with /if-band magnitude increases to 97% 
in Wl and 95% in W4. 

We compare the K -band photometry for optical sources 
matched with both UKIDSS and WIRCAM surveys, and find 
a good agreement. In fact, we find a mean difference (AK) = 
(^wircam - #ukidss> - -0.05, with a small dispersion cr AK 
0.10 and 0.15, for Wl and W4 respectively. These differences 
can be ascribed to the transmission functions of the filters and 
the definition of the aperture used when measuring magnitudes, 
and are of the order of photometric errors. To not overweight the 
K-bmd magnitudes in the SED fitting, only the deeper /ifwiRCAM 
has been used when both magnitudes were available for the same 
object. 



3 DR9 for LAS and DXS, DR8 for UDS; http : //www . ukidss . org/ 



Besides the NIR magnitudes, the UV part of the spectrum 
can also be important for constraining the galaxy dust content 
and the star formation rate. We make use of existing GALEX 
images observed with the deep imaging survey (integration time 
~ 3 x 10 4 s) in NUV and FUV channel and we have completed 
the coverage in Wl region with new observations in NUV chan- 
nel only and with integration time T exp > 1.5 x 10 4 s. Because 
of the GALEX large PSF (~ 5 arcsec), the source blending is 
a major issue in GALEX deep imaging mode. To measure the 
UV fluxes of the sources, we use the dedicated photometric algo- 
rithm EMphot (Conseil et al. 201 1) which adopts the positions of 
U band selected priors and performs a modelled PSF adjustment 
over small tiles based on the expectation maximisation algorithm 
(Guillaume et al. 2006). For our spectroscopic sample, 63% 
(15%) of the sources have a NUV (FUV) flux measurements 
in Wl. In contrast, the W4 field has a modest GALEX coverage 
with 13% (5%) of spectroscopic sources with a NUV (FUV) 
flux. The WIRCAM and GALEX datasets in the VIPERS fields 
are described in Arnouts et al. (in prep.). 

Moreover, for ~ 30% of the spectroscopic targets already 
observed in Wl, we also took advantage of the SWIRE obser- 
vations in the XMM-LSS field. During the SED fitting we only 
considered magnitudes in the 3.6 //m and 4.5 fim bands, since at 
longer wavelengths the source detection is very sparse. More- 
over, at those wavelengths the re-emission from the dust begins 
to contribute to the flux of galaxies, and this feature is not repro- 
duced by most of the models of stellar population synthesis (see 
Sect. 2.3). 

2.2. Spectroscopy 

The spectroscopic catalogue used in this paper represents the 
first 60% of VIPERS. This sample includes 53608 galaxy spec- 
tra and will be made available through the future VIPERS Pub- 
lic Data Release 1 (PDR-1). The VIPERS targets were selected 
via two criteria. The first was aimed at separating galaxies and 
stars, and relies on the combination of a point-like classification 
(based on the measure of the half-light radius) for the bright- 
est sources, and on the comparison of the five optical magni- 
tudes with galaxy and stellar spectral energy distributions for 
the faintest ones (Coupon et al. 2009). A fraction of the point- 
like sources are targeted as AGN candidates, when located in the 
AGN loci of the two colour diagrams (g - r) versus (u - g) and 
(g - i) versus (u - g). The second selection criterion, based on 
(g - r) and (r - i) colours, has been developed to exclude low 
redshift (z < 0.5) objects, and has been tested to ensure it does 
not introduce any significant bias. A complete description of 
the whole source selection procedure is included in Guzzo et al. 
(2013). 

The spectroscopic observations have been carried out using 
the VIMOS instrument on VLT LR-Red grism (R = 210), giv- 
ing a wavelength range of 5500-9500 A that guarantees the ob- 
servability of the main spectral features in the VIPERS redshift 
range, e.g. the absorption lines Call H & K /L13934, 3969 and 
the emission line [Oil] A3721 . Using a sample of objects spec- 
troscopically observed twice, we are able to estimate an uncer- 
tainty of cr z = 0.00047 (1 + z) for our measured redshifts. 

To maximise the multiplex capability of VIMOS, we adopt 
the observational strategy described in Scodeggio et al. (2009) of 
using shorter slits than in the previous surveys carried out with 
the same instrument. By virtue of this strategy, we reached a 
sampling rate of approximately 40% with a single pass, essential 
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to estimate the large scale environment (Cucciati et al., in prep.; 
Iovino et al., in prep). 

The spectroscopic masks reproduce the footprint of the VI- 
MOS instrument, consisting of 4 quadrants and gaps between 
them for each pointing, covering 224arcmin 2 . Moreover, vi- 
gnetted parts of the quadrants have been removed to compute the 
effective area (for a detailed description see Guzzo et al. 2013). 

Data reduction and redshift measurement have been per- 
formed within the software environment Easylife (Garilli et al. 
2012), based on the VIPGI pipeline (Scodeggio et al. 2005) and 
EZ (Garilli et al. 2010, Easy redshift). Spectroscopic redshifts 
have been determined independently by two team members and 
reconciled in case of discordance. Each redshift has a confidence 
assigned by using a well-established schemed developed by pre- 
vious surveys like VVDS (Le Fevre et al. 2005) and zCOSMOS 
(Lilly et al. 2009). In detail, a spectroscopic quality flag of 4 
corresponds to a confidence level of 99.6% (Guzzo et al. 2013), 
down to a confidence level of 0, which represents the lack of a 
reliable redshift estimate. Moreover, objects with a single emis- 
sion line are labelled by flag 9, and broad-line AGNs share the 
same scheme, but their flags are increased by 10. Each spec- 
troscopic flag also has a decimal digit specifying the agreement 
with the photometric redshift computed from CFHTLS photom- 
etry (Coupon et al. 2009). 

After excluding galaxies without redshifts (flag 0) and stars, 
our redshift sample contains 53 608 extragalactic sources, nearly 
equally split between the two fields. The quality of redshift mea- 
surements for the sample with spectroscopic flags larger than 2, 
as estimated from the validation of multiple observations, is high 
(confidence > 95%, see Guzzo et al. 2013). 

Since a fraction of all the possible targets have been ob- 
served, statistical weights are required to make this subsample 
representative of all the galaxies at i < 22.5 in the survey vol- 
ume. Such weights are calculated taking into account the num- 
ber of photometric objects that have been targeted (target sam- 
pling rate, TSR), the fraction of them classified as secure mea- 
surements (spectroscopic success rate, SSR), and the complete- 
ness due to the colour selection (colour sampling rate, CSR; for 
a detailed discussion see Guzzo et al. 2013; Fritz et al. 2013). 
The statistical weights can depend on the magnitude, redshift, 
colour, and angular position of the considered object. For each 
part of the statistical weight we considered only the main and 
relevant dependencies, in order to avoid spurious fluctuations in 
case of small subsamples. 

In particular, we considered the TSR as a function only of 
the selection magnitude, the SSR as a function of magnitude and 
redshift, and the CSR (estimated by using data from the VVDS 
flux limited survey, Le Fevre et al. 2005) as a function of redshift. 
Regarding the SSR, only galaxies with quality flag between 2 
and9(~ 41 100 galaxies in the redshift range 0.5 < z =S 1.3) have 
been considered in the analysis (we exclude spectra classified as 
broad-line AGNs). For a galaxy at redshift z with magnitude its 
statistical weight w(i, z) is the inverse of the product of TSR(0, 
SSR(/, z), and CSR(z). Once each galaxy in the spectroscopic 
sample is properly weighted, we can recover the properties of the 
photometric parent sample with good precision (see Fritz et al. 
2013, for the details on TSR, SSR, and CSR). 

2.3. Stellar masses 

Considering the small fraction of objects without K band mag- 
nitude, we decided to rely on SED fitting to derive stellar masses 
and to not implement alternative methods such as the Lin et al. 
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Fig. 1. Distribution of the differences between the values of stel- 
lar mass estimated using the two codes Hyperzmass and MAGPHYS. 
Only results on the Wl field are show (see text). To limit the effects 
of parameter degeneracy, we restrict the comparison to galaxies that 
turn out to have solar metallicity, according to their best- fit SED both in 
Hyperzmass and MAGPHYS. In this way log(WI M AGPHYs/AlHyperzma SS ) 
cannot be due to a different metal content assumed in the two SED 
fitting estimates. The red dashed line gives the best-fitting Gaussian 
of the distribution, corresponding to the mean and standard deviation 
indicated. Also indicated are the size of the galaxy subsample (N m ) 
and the number of stellar mass estimates for which the discrepancy is 

log( MmAGPH YS / ^Hypei ) > 2(T (iVexceed)- 

(2007) relation between stellar mass, redshift and rest-frame 
magnitudes. 

We thus derive galaxy stellar masses ( At) by means of an up- 
dated version of Hyperzmass (Bolzonella et al. 2000, 2010, the 
software is available on request). Given a set of synthetic spec- 
tral energy distributions, the software fits these models to the 
multi-band photometry for each galaxy and selects the model 
that minimises the^f 2 . The SED templates involved in this pro- 
cedure are derived from simple stellar populations (SSP) mod- 
elled by Bruzual & Chariot (Bruzual & Chariot 2003, hereafter 
BC03), adopting the Chabrier (2003) universal initial mass func- 
tion (IMF) 4 . The BC03 model is one of the most used in lit- 
erature (e.g. Ilbert et al. 2010; Zahid et al. 2011; Barro et al. 
2013). Another commonly used SSP library is the one by Maras- 
ton (M05, 2005), which differs from the former because of the 
treatment of the thermally pulsing asymptotic giant branch (TP- 
AGB) stellar phase, especially affecting NIR emission of stellar 
populations aged ~ 1 Gyr. The issue about the TP-AGB role 
in the stellar population synthesis is still open (e.g. Marigo & 
Girardi 2007), with some evidences that supports BC03 (Kriek 
et al. 2010; Zibetti et al. 2013) in contrast with observations 
favouring M05 (MacArthur et al. 2010). In the following we 
prefer to adopt the BC03 model, since most of the galaxies 
in the redshift range we explored should not be dominated by 
the TP-AGB phase (which is instead relevant for galaxies at 
1.4 < z < 2.7 Maraston et al. 2006). 



4 The choice of a different IMF turns into a systematic mean offset 
in the stellar mass distribution: for instance, our estimates can be con- 
verted to a Salpeter (1955) IMF by a scaling factor of ~ 1.7, and ~ 1.1 
for switching to Kroupa (2001). 
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The SSPs provided by Bruzual & Chariot (2003) assume a 
non-evolving stellar metallicity Z, which we chose to be solar 
(Z = Z G ) or sub-solar (Z = 0.2 Z Q ). This choice allows us to 
take into account the different metallicities of the galaxies in our 
redshift range, which can be lower than in the nearby universe 
(Zahid et al. 2011), without significantly increasing the effect 
of the age-metallicity degeneracy. Considering the low resolu- 
tion of our spectroscopic setup, it is difficult to put reliable con- 
straints on Z from the observed spectral features, and therefore it 
was not possible to constrain this parameter a priori. Therefore, 
the metallicity assigned to each galaxy is that obtained from the 
best-fit model (smallest^ 2 ). 

With respect to the galaxy dust content, we implement the 
Calzetti et al. (2000) and Prevot-Bouchet (Prevot et al. 1984; 
Bouchet et al. 1985) extinction model, with values of Ay rang- 
ing from (no dust) to 3 magnitudes. As pointed out by Ilbert 
et al. (2009), Calzetti's law is on average more suitable for the 
bluest SEDs, having been calibrated on starburst (SB) galaxies, 
whereas the Prevot-Bouchet law is better for mild star-forming 
galaxies, as it was derived from the dust attenuation of the Small 
Magellanic Cloud (SMC) (see also Wuyts et al. 201 1). Hereafter 
we refer to the Calzetti and Prevot-Bouchet models as SB and 
SMC extinction laws respectively. We let the choice between the 
two extinction laws free, according to the best fit model (small- 
est x 2 ), since we do not have sufficient data at UV wavelengths 
to discriminate between the different trends of the two laws. 

The SEDs constituting our template library are then gener- 
ated from the SSPs following the evolution described by a given 
star formation history (SFH). In this work, we assume expo- 
nentially declining SFHs, for which SFR <x exp(-?/r), with the 
time scale t ranging from 0.1 to 30Gyr. A constant SFH (i.e., 
SFR ~ 1 Aloyr -1 ) is also considered. 

Although such a parametrisation is widely used, recent stud- 
ies have shown how exponentially increasing SFHs can provide a 
more realistic model for actively star-forming galaxies in which 
young stellar populations out-shine the older ones (Maraston 
et al. 2010). This effect becomes relevant at z ~ 2, when the 
cosmic star-formation peaks, and can be reduced by setting a 
lower limit on the age parameter, in order to avoid unrealistic 
solutions that are too young and too dusty (Pforr et al. 2012). In 
our redshift range, galaxies whose SFH rises progressively have 
small stellar masses (\og(M/M B ) ~ 9.5, Pacifici et al. 2013) 
falling below the limit of VIPERS. 

Moreover, Pacifici et al. (2013) identify a class of massive 
blue galaxies that assembled their stellar mass over a relatively 
long period, experiencing at a later stage the progressive reduc- 
tion of their star formation. For such bell-shaped SFH, neither 
increasing nor decreasing r-models seem to be suitable. How- 
ever, the resulting differences are smaller than the other uncer- 
tainties of the SED fitting method (cf. Conroy et al. 2009), and 
the stellar masses computed assuming r-models are adequate for 
our analysis. 

Another issue concerning the SFH is the assumption of its 
smoothness. In fact, a real galaxy could have experienced sev- 
eral phases of intense star-formation during its past, which can 
be taken into account by superimposing random peaks on the ex- 
ponential (or constant) SFR (Kauffmann et al. 2003a). Allowing 
the presence of recent secondary bursts, making the colours of 
an underlying older and redder population bluer, can lead to a 
systematically larger stellar mass estimate, but only for a small 
fraction of objects the difference in M is larger than 0.2 dex, as 
shown by Pozzetti et al. (2007). 

We also quantified the effect of complex SFHs in VIPERS, 
by computing stellar masses using the MAGPHYS package (da 



Cunha et al. 2008). This code parametrises the star formation ac- 
tivity of each galaxy template starting from the same SSP mod- 
els as Hyperzmass (i.e. BC03), but using two components in the 
SFH, namely an exponentially declining SFR and a second com- 
ponent of additional bursts randomly superimposed to the former 
according to Kauffmann et al. (2003a). The probability of a sec- 
ondary burst to occurring is such that half of the galaxy templates 
in the library have experienced a burst in their last 2 Gyr. Each 
of those episodes can last 3 x 10 7 — 3 x 10 8 yr, producing stars 
at a constant rate. The ratio between the stellar mass produced 
in a single burst and that formed over the entire galaxy life by 
the underlying exponentially declining model is distributed log- 
arithmically between 0.03 and 4.0. MAGPHYS also treats dust 
attenuation in a consistent way, adopting the model proposed by 
Chariot & Fall (2000), and with the re-emission at infrared wave- 
lengths. However, this feature does not represent a significant 
advantage when dealing with VIPERS data as infrared magni- 
tudes are too sparse in our catalogue. Metallicity values are dis- 
tributed uniformly between 0.02 and 2Z . The large range of 
tightly sampled metallicities, the different model for the dust ex- 
tinction, and in particular the complex SFHs in the MAGPHYS 
library are the major differences to the Hyperzmass code. 

In Fig. 1 we compare the estimates obtained through MAG- 
PHYS and Hyperzmass, and verify that complex SFHs have 
a minimal impact on the results. Since MAGPHYS requires 
much longer computational time than other SED fitting codes, 
we only estimate the stellar mass for galaxies in the Wl field 
between z = 0.5 and z = 1.3. Moreover, we select objects 
with the same (solar) metallicity in both the SED fitting pro- 
cedures. In this way we are able to investigate the bias mainly 
due to the different SFH parametrisations. We find a small off- 
set ((AlogM) = <log(AWcPHYs/AlHypei-zmas S )> - 0.05) and a 
small dispersion (<t&m — 0.11), with significant differences be- 
tween MAGPHYS and Hyperzmass (i.e., AlogAI > 0.22) for 
~ 7% of the testing sample (A exC eed in Fig. 1). The consequences 
on the GSMF are discussed in Sect. 3.4. 

Given the large range of physical properties allowed in the 
SED fitting procedure, we decide to exclude from the fitting 
some unphysical parameter combinations. In particular, we limit 
the amount of dust in passive galaxies (i.e., we impose Ay < 0.6 
for galaxies with age/r > 4), we avoided very young extremely 
star-forming galaxies with short t timescales (i.e. we prevent 
fits with models with r < 0.6 Gyr when requiring Zf 0lm < 1), and 
we only allow ages to be within the range between 0.1 Gyr and 
the age of the universe at the spectroscopic redshift of the fitted 
galaxy (see Pozzetti et al. 2007; Bolzonella et al. 2010). 

According to Conroy et al. (2009), the uncertainties associ- 
ated with the SED fitting can be ~ 0.3 dex when considering all 
the possible parameters involved and their allowed ranges. In 
particular, given the non-uniform coverage of the GALEX and 
SWIRE ancillary data matched with our sample, we checked that 
the variation in the magnitude set from one object to another is 
not introducing any bias. For the subsample of galaxies with 
FUV, NUV, 3.6pm, and 4.5 pm bands available, we estimate 
the stellar mass also using just the the optical-NIR photometry. 
We find no systematic difference in the two estimates of stellar 
mass (with and without the UV and infrared photometry) and 
only a small dispersion of about 0.08 dex. 

In summary, the VIPERS galaxy stellar mass estimates 
are obtained using the BC03 population synthesis models with 
Chabrier IMF, smooth (exponentially declining or constant) 
SFHs, solar and sub-solar metallicity, and the SB and SMC laws 
for modelling dust extinction. Unless stated otherwise, this is the 
default parametrisation used throughout this paper. 
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3. From stellar masses to the galaxy stellar mass 
function 

In this Section we exploit the VIPERS dataset described above, 
considering only our fiducial sample of 41 094 galaxies at z — 
[0.5,1.3] with spectroscopic redshift reliability > 95% (see 
Sect. 2.2). As mentioned before, broad line AGN (~ 850 in the 
present spectroscopic sample) are naturally excluded from the 
sample, being visually identified during the redshift measure- 
ment process. Instead, narrow line AGN are not removed from 
our sample, but they do not constitute an issue for the SED fit- 
ting derived properties, since in most of the cases their optical 
and NIR emission are dominated by the host galaxy (Pozzi et al. 
2007). First of all we try to identify the threshold below which 
the sample can be considered complete, and the mass function 
reliable. After that, we derive the GSMF of VIPERS in various 
redshift bins and discuss the main sources of uncertainty affect- 
ing it. 

3.1. Completeness 

In the literature, the completeness mass limit of a sample at a 
given redshift is often defined as the highest stellar mass a galaxy 
could have, when its observed magnitude matches the flux limit 
(e.g. Perez-Gonzalez et al. 2008). In that case, the maximum of 
the stellar mass estimates is reached by the rescaled SED of an 
old passive galaxy, which gives rise to a limiting mass that is 
too conservative. In fact, because of the variety of mass-to-light 
ratios Ai/L characterising the different galaxy types, this bound- 
ary in stellar mass strongly depends on the observed galaxy 
mix. At the limiting magnitude of a survey, blue galaxies (with 
low Ai/L) are the dominant population (e.g. Zucca et al. 2006), 
whereas red galaxies are generally less numerous. Therefore, 
the Ai/L range of galaxies around the magnitude limit is much 
narrower than the same range for the total sample, and usually 
it does not include the large Ai/L values (typical of old passive 
galaxies) used by the authors to estimate the completeness (see 
the discussion in Marchesini et al. 2009, Appendix C). 

Instead, we apply the technique devised by Pozzetti et al. 
(2010). This procedure yields, at a certain redshift and for a 
fixed flux limit, an estimate of the threshold Ain m below which 
some galaxy type cannot be detected any longer. Following this 
approach, we estimate the stellar mass each object would have 
if its magnitude, at the observed redshift, were equal to the i- 
band limiting magnitude ;'ij m . This boundary mass Ai(i = him) is 
obtained by rescaling the original stellar mass of the source at its 
redshift, i.e. log Ai(i = in m ) = log Ai + 0A(i - iu m ). The threshold 
Alum is then defined as the value above which 90% of the Ai(i = 
iiim) distribution lies. According to this, at values larger than 
Alum our GSMF can be considered complete. We include in 
the computation only the 20% faintest objects to mitigate the 
contribution of the small number of bright red galaxies having 
large Ai/L. In this way we improve the method mentioned at the 
beginning of this Section, allowing us to extend the measurement 
to lower masses. 

Since the 1/V ma x method (Schmidt 1968, see Sect. 3.2) in- 
trinsically corrects the sample incompleteness above the lower 
limit of the considered redshift bin Zwi , we apply to each redshift 
bin the Aiy, m computed considering the objects inside a narrow 
redshift interval Az = 0.05 centred in Zj n f . Figure 2 shows A1ii m 
as a function of redshift for the global and for the red and blue 
samples used in Sect. 5, as well as the value of Ai(i = ihm) for 
each red and blue galaxy. As expected, the limiting mass in- 
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Fig. 2. The mass completeness threshold Mn m as a function of red- 
shift, computed for the total sample (the one used in Sect. 3.2, filled 
circles) and for the red (upward triangles) and blue (downward trian- 
gles) populations, defined as discussed in Sect. 5. In each redshift bin, 
the M\\ m estimate relies on the rescaled stellar mass Mil = i]\ m ) of the 
20% faintest galaxies (see text). We show M(i = iW) of the red and 
blue galaxies with small dots of analogous colours. 
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Fig. 3. Distributions of stellar masses in six redshift bins for the 
VVDS-Deep sample in the CFHTLS-W1 field at its limiting magnitude 
(/ < 24, dark histograms), compared to the subset obtained by applying 
a VIPERS-like magnitude limit / < 22.5 (blue histograms). In each 
panel, the black dashed line represents the limiting mass for the latter 
sample, according to the method described in Sect. 3.1. The red solid 
line instead gives the limiting mass estimated through the same method, 
but for the VIPERS sample in the Wl field. Both limits correctly iden- 
tify the threshold below which the shallower sample starts to miss a 
significant fraction (> 20%) of objects. 



creases as a function of z and the values for red galaxies are 
much larger than for the blue ones. 

In the context of the zCOSMOS project (Lilly et al. 2009), 
the approach of Pozzetti et al. (2010) produced completeness 
limits in good agreement with those obtained through mock sur- 
vey samples (Meneux et al. 2009). In VIPERS, we successfully 
tested our Ai\i m estimates taking advantage of the VVDS-Deep 
field, which is located in the Wl field (see Guzzo et al. 2013, 
Fig. 2). The VVDS sample provides us with spectroscopically 
observed galaxies up to a fainter limit, /ab = 24 (Le Fevre et al. 
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Fig. 4. The VIPERS galaxy stellar mass function at different redshifts. The circles give the values determined through 1/V max in mass bins of 
AM = 0.2 dex; the centre of each bin corresponds to the weighted mean mass of the objects within it. Empty and filled symbols correspond to 
values below and above the completeness limit, respectively (see Sect. 3.1). For the latter points, the red error bars show the uncertainty due to 
Poisson noise, while green bars account for Poisson noise and cosmic variance. In each panel, a solid line shows the Schechter fit to the GSMF 
points, with the dashed line reproducing that of the first redshift bin, as a reference. The arrows give an upper limit to C> where no detection is 
available. 



2005). Since the CFHTLS-W1 field contains both VVDS and 
part of VIPERS, we can compare the stellar masses relying on 
a similar photometric baseline {u, g, r, I, i, z, J*, K*). When ap- 
plying a VIPERS-like magnitude cut (/ < 22.5), we can find the 
fraction of missed objects with respect to the parent / < 24 sam- 
ple as a function of stellar mass. This test is shown in Fig. 3, 
where we compare the values of the VVDS and VIPERS M\\ m 
for / < 22.5 to the distribution of stellar masses belonging to 
the deeper VVDS sample. The Wliim values we computed rep- 
resent well the limit at which the stellar mass distribution starts 
to be incomplete with respect to the deep VVDS sample (i.e. the 
/ < 22.5 sample recovers less than 80% of the parent sample). 

3.2. Evolution of the mass function for the global population 

The number of galaxies and volume sampled by VIPERS allows 
us to obtain an estimate of the GSMF with high statistical preci- 
sion within six redshifts bins in the range 0.5 < z < 1.3. Given 
the large number of galaxies observed by VIPERS, in terms of 
Poisson noise it would be possible to choose even narrower bins 
(e.g. Az - 0.05 wide). However, in that case the measurements 
start being strongly affected by cosmic (sample) variance. A 
more detailed discussion is given in Sect. 3.3. 

We compute the GSMF within each redshift bin, using the 
classical non-parametric 1 /V max estimator (Schmidt 1968). With 
this method, the density of galaxies in a given stellar mass bin is 
obtained as the sum of the inverse of the volumes in which each 
galaxy would be observable, multiplied by the statistical weight 
described in Sect. 2.2. To optimise the binning in stellar mass, 



we use an adaptive algorithm that extends the width of a bin until 
it contains a minimum number of objects. The errors associated 
with the 1 /V max estimates are computed assuming Poisson statis- 
tics. The upper limits for non-detections have been estimated 
following Gehrels (1986). The values of the 1/V max GSMF and 
associated Poisson errors are given in Table 1 . 

It is well known that the 1 / V max estimator is unbiased in case 
of a homogeneous distribution of sources (Felten 1976), but it 
is affected by the presence of clustering (Takeuchi et al. 2000). 
At variance with the data sets on which the estimator was tested 
in the past, VIPERS has a specific advantage, thanks to its large 
volume over two independent fields. The competing effects of 
over- and under-dense regions on the estimate should in such 
a situation cancel out. Moreover, as shown in Takeuchi et al. 
(2000), an inhomogeneous distribution of sources mainly affects 
the faint end (i.e. the low mass end) of the luminosity (mass) 
function, and the impact on our analysis will be negligible, as 
we are mainly interested in the massive tail of the distribution. 

To verify this, we also compare the 1/V max estimates with 
those from a different software package (ALF, Ilbert et al. 2005), 
which includes several estimators for the luminosity and mass 
functions. We find no significant differences in the obtained 
mass functions, within the stellar mass range considered in the 
present study. 

Finally, in addition to the non-parametric method, we fit a 
Schechter function (Schechter 1976) to the 1 /V max estimates: 
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Table 1. Global GSMF: 1/V max values in regular bins of stellar mass. 
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The results are shown in Fig. 4. The 1 / V max determinations 
are marked with filled circles when they are above the complete- 
ness limit Aliim (see Sect. 3.1). Although the mass function does 
not show evidence for a rapid decline below this threshold (cf. 
Drory et al. 2009), the open circles should be conservatively con- 
sidered as lower limits. The solid lines give the best-fit using a 
Schechter function. 

These plots show clearly the statistical power of the VIPERS 
sample, which includes a significant number of the rare massive 
galaxies that populate the GSMF high-mass end, thanks to its 
large volume. 

At z < 0.6 we detect the characteristic dip of the mass func- 
tion at log(A1/A1o) ~ 10.2, with an upturn below that value as 
observed both locally (e.g. Baldry et al. 2012) and at intermedi- 
ate redshifts (e.g Drory et al. 2009; Pozzetti et al. 2010). How- 
ever, we avoid using a double Schechter function in our fits, to 
ease comparison with the parameters derived at higher redshifts. 
In fitting the points in the first bin (0.5 < z < 0.6), all parameters 
of Eq. 1 are left free, obtaining a value of the slope a = -0.95. 
Above this redshift, however, the slope of the low-mass end is 
only weakly constrained, given the relatively large values of the 
completeness limit Alii m . For this reason, in all the other bins 
we fix a to the value -0.95 (see Table 2). 

The results of Fig. 4 show with impressive statistical pre- 
cision the lack of evolution since z - 1.1 of the massive end 
(log(A1/A1o) > 11) of the galaxy mass function seen in previ- 
ous, smaller samples. The exponential tail of the Schechter fit 
is in fact nearly constant across the five redshift bins, down to 
z — 0.5 (see Fig. 4). However, we detect a significant decrease 
between z = 1.1 and 1.3, of the number density over the same 
mass range. At lower masses (10.8 < log(A1/A1o) < 11.1) the 
first signs of evolution with respect to z ~ 0.5 start to be visible 
at redshift 0.9- 1.1. 

These trends can be better evidenced in Fig. 5, where the 
number density of galaxies p# within three mass ranges is plot- 
ted versus redshift. This figure explicitly shows that the most 
massive galaxies are virtually already in place at z - 1 . On the 
contrary, galaxies with smaller mass keep assembling their stars 
in such way that their number density increases by a factor ~ 3.5 



Table 2. Global GSMF: Schechter parameters 
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from z = 1.2 down to 0.6, consistently with the so-called down- 
sizing scenario (Cowie et al. 1996; Fontanot et al. 2009). These 
new measurements confirm previous evidence, but with a higher 
statistical reliability. In particular, as we shall discuss in Sect. 4, 
these estimates have a precision, in terms of Poisson error, that 
was never reached before at such redshifts by a spectroscopic 
survey. This improvement is paralleled also by other sources of 
error, in particular cosmic variance, which we shall discuss in 
detail in Sects. 3.3 and 3.4. 



3.3. Cosmic variance in the VIPERS survey 

When dealing with statistical studies using number counts, a se- 
vere complication is introduced by the field-to-field fluctuations 
in the source density, due to the clustered nature of the galaxy 
distribution and the existence of fluctuations on scales compa- 
rable to the survey volume. This sampling or 'cosmic' variance 
represents a further term of uncertainty to be added to the Pois- 
son shot noise. It can be expressed by removing cr 2 oiss si/ (N) 
from the total relative error: 



(jV 2 ) - (N) 2 



1 



(2) 
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Fig. 5. Evolution of the galaxy number density in different bins of stel- 
lar mass. The error bars of the density estimates includes Poisson noise 
and cosmic variance (see Sect. 3.3). At z - 1.2, for the lowest-mass 
sample, only a lower limit can be estimated, indicated by the arrow. 

where (N) and (A^ 2 ) are the mean and the variance of galaxy 
number counts (Somerville et al. 2004). 

Extragalactic pencil-beam surveys, even the deepest ones, 
are particularly limited by cosmic variance, given the small vol- 
ume covered per redshift interval. At z ~ 0.8, galaxy den- 
sity fluctuations are found to be still relevant up to a scale of 
~ 140Mpc/i^ ( J (Scrimgeour et al. 2012), which roughly corre- 
sponds to 5 deg. 

This is the result of intrinsic clustering in the matter, as pre- 
dicted by the power spectrum shape and amplitude at that epoch, 
amplified by the bias factor of the class of galaxies analysed, 
which at high redshift can be very large for some classes. For 
these reasons, even the largest deep surveys of the past genera- 
tion do not allow us to overcome this issue. For example, the 
COSMOS field, despite its 2 deg" area, turned out to be signifi- 
cantly over-dense between z = 0.8 and z = 1 (Kovac et al. 2010). 

The gain obtained by enlarging the area of a single field 
beyond a certain coverage becomes less prominent, due to the 
existing large-scale correlations (see Newman & Davis 2002, 
Fig. 1): <x cv decreases mildly as a function of volume, with an 
approximate dependence cr cv oc y~° 3 (Somerville et al. 2004, 
Fig. 2), compared to cr Po j ss oc y~ - 5 . Trenti & Stiavelli (2008) 
found similar results by characterizing Lyman break galaxies 
surveys: at high values of (AO, the Poisson noise rapidly drops 
and cosmic variance remains the dominant source of uncertainty. 
A more effective way to abate cosmic variance is to observe sep- 
arated regions of sky. Since counts in these regions, if they are 
sufficiently distant, are uncorrected, their variances sum up in 
quadrature (i.e., cr cv decreases as the square root of the number 
of fields, Moster et al. 2011). Multiple independent fields can 
then result in an uncertainty smaller than a single one, even if 
the latter has a larger effective area (Trenti & Stiavelli 2008). 
The current VIPERS PDR-1 sample is not only characterized by 
a significantly large area (thus volume), compared to previous 
similar surveys at these redshifts, but this is also split into two 
independent and well-separated fields of ~ 7.5 deg 2 each. We 
therefore expect that the impact of cosmic variance should be 
limited. 



To quantify the impact of cosmic variance directly, we fol- 
low two approaches. The first one, based on the observations 
themselves, provides an upper limit of the VIPERS cr cv . We 
select five rectangular subregions of about 2 deg 2 within the sur- 
vey and estimate the mass function O, in each of them, using the 
1 / Vmax method described above. We choose non-contiguous re- 
gions (separated by ~ 1 deg) to minimize the covariance between 
subsamples located within the same field (Wl or W4). Within 
mass bins Mj + AA1/2 we derive the total random uncertainty 



cr tot , obs (Mj) = i £ - <t> tot (Mj)f , (3) 

1=1 

where <D tot is the global GSMF of VIPERS (at that redshift) and 
<t>,(AL/) is the number density of galaxies measured in the j-th 
mass bin for each of the n = 5 subregions. 

The second approach is based on the use of simulated mock 
surveys. First, we use a set of 57 mock samples (26 and 3 1 in Wl 
and W4, respectively), built using specific recipes for the stellar- 
to-halo mass relation. They are based on the MultiDark dark 
matter simulation (Prada et al. 2012) and have been constructed 
to reproduce the detailed geometry and selection function of the 
VIPERS survey (see de la Torre et al. 2013, for details). In that 
process, the dark matter haloes identified in the simulation, as 
well as artificial sub-haloes drawn from the Giocoli et al. (2010) 
sub-halo mass function, have been associated with galaxies us- 
ing the stellar-to-halo mass relations of Moster et al. (2013). The 
latter are calibrated on previous stellar mass function measure- 
ments in the redshift range < z < 4. We call these 'SHMR 
mocks'. We apply Eq. 2 to estimate the amount of cosmic vari- 
ance independently among the 26 Wl and 31 W4 mocks. The 
global estimate of cosmic variance (<x cv shmr) on the scales of 
the VIPERS survey is obtained by combining the results from 
the two fields (see Moster et al. 2011, Eq. 7). As expected, we 
find that cr cv SHMR decreases with redshift, since we are probing 
larger and larger volumes, and increases with stellar mass due 
to the higher bias factor (and thus higher clustering) of massive 
galaxies (Somerville et al. 2004). Both the trends are clearly vis- 
ible in Fig. 6 , where measurements of <t cv shmr are presented 
for different bins of redshift and stellar mass. 

Besides these SHMR mocks, we also used another set of 50 
VIPERS-like light-cones built from the Millennium simulation 
(Springel et al. 2005), in which dark-matter haloes are populated 
with galaxies through the semi-analytical model (SAM) of De 
Lucia & Blaizot (2007). Galaxy properties were determined by 
connecting the astrophysical processes with the mass accretion 
history of the simulated dark matter haloes. Each mock sample 
covers 7x1 deg 2 , with a magnitude cut in the i band equal to 
that of the observed sample. Although the geometry of these 
mocks (and therefore their volume) slightly differs from the de- 
sign of the real survey, they provide an independent test, with 
completely different prescription for galaxy formation, with re- 
spect to the SHMR mocks. 

In Fig. 6 (crosses), we can notice that the estimates of the 
standard deviation obtained from the field-to-field fluctuations 
among the five sub-samples, cr tot obs have an erratic behaviour 
and do not show the smooth trend indicated by of cr cv SHMR . The 
squares should in fact be regarded as an upper limit, given that 
the subsamples have a volume smaller than the whole survey and 
since Eq. 3 includes also the variance due to Poisson noise. Con- 
versely, residual correlation among the sub-fields within each of 
the VIPERS fields (produced by structures on scales > 1 deg 
crossing over two or more subregions) would reduce cr tot obs . 
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Fig. 6. Estimates of the contribution of sample (cosmic) variance 
to the statistical uncertainty of the GSMF measurements. For each 
redshift bin, the upper panels show the GSMF 1/V m ax measurements 
obtained from five 2 deg 2 VIPERS subregions, located respectively in 
the Wl field (three regions, blue diamonds, circles and squares) and in 
the W4 field (two regions, red triangles and downward triangles). The 
Schechter fit to the global GSMF of Fig. 4 is shown as reference (black 
solid line). The lower panels report the standard deviations estimated in 
each redshift bin from these five measurements (purple squares, Eq. 3), 
together with the estimates of cr cv obtained from 57 SHMR mocks by 
means of Eq. 2 (green crosses). To highlight how the effect of cosmic 
variance decreases at higher z, we report o" cv SH mr of the first redshift 
bin in the other panels. In addition, the sample variance measured in 
50 SAM mocks (grey solid line) and the estimates provided by Moster 
et al. (201 1) method (black triangles) are shown as reference. 



More in general, the small number of fields used to perform this 
test makes the computation of Eq. 3 statistically uncertain. 

In that figure, we also show, as a reference, the estimates pro- 
vided by the public code getcv (Moster et al. 201 1) for the same 
area and redshift bins of the SHMR mocks; these results, limited 
at log(A1/Al Q ) < 1 1 .5, are in good agreement with cr cv SHMR . 



3.4. Other sources of uncertainty 

In describing our procedure to derive stellar masses by means of 
the SED fitting technique (Sect. 2.3), we emphasised the number 
of involved parameters, and their possible influence on the esti- 
mates. The assumptions that have the strongest impact on the 
results are the choices of the stellar population synthesis model, 
IMF, SFH, metallicity, and dust extinction law. A thorough dis- 
cussion about each one of the mentioned ingredients is beyond 
the goals of this paper, and the reader can refer to Conroy (2013) 
and Marchesini et al. (2009) for a comprehensive review of the 
systematic effects induced by the choice of the input parameters. 

Here we briefly test the impact on the GSMF of choosing 
different values of Z (whether including sub-solar metallicities 
or not), the extinction laws (SB and SMC, or SB alone), and the 
addition of secondary bursts to the smooth SFHs (i.e. complex 
SFHs instead of exponentially declining r-models). We will not 
modify the other two main ingredients in our procedure, i.e. the 
universal IMF we assumed (Chabrier 2003) and the stellar pop- 
ulation synthesis model (BC03). 

To perform this test we use stellar mass estimates obtained 
by assuming five different sets of SED fitting templates, four of 
them differing in metallicity and extinction law: Z only and 
SB; two metallicities (Z and O.2Z ) and SB; solar metallicity 
and two extinction laws (SB and SMC); two metallicities (Z 
and O.2Z ) and two extinction laws (SB and SMC). The fifth 
SED fitting estimate has been derived with the MAGPHYS code 
(see Sect. 2.3), assuming the following parameters: complex 
SFHs, extinction model derived from Chariot & Fall (2000), and 
a larger range of metallicity (including super-solar ones) . We 
limit these tests to the data in the VIPERS Wl field, i.e. about 
half of the total sample, given the better overall photometric cov- 
erage in this area and the large computational time involved. 

The mass functions resulting from these five different SED- 
modelling assumptions are reported in Fig. 7. As expected (see 
discussion in Sect. 2.3 and Fig. 1), the MAGPHYS mass function 
corresponds to the highest estimated values of galaxy density 
at large stellar masses (at least up to z - 1.1). This trend is not 
unexpected, because the other four estimates, obtained assuming 
smooth SFHs templates, are insensitive to an underlying old stel- 
lar population that is outshone by a recent burst of star formation 
(Fontanaet al. 2004; Pozzetti et al. 2010, but see Moustakas et al. 
2013 for an opposite result). As a consequence, using complex 
SFHs templates one can produce stellar mass estimates that are 
larger than those obtained with smooth SFHs for a small percent- 
age of objects, an effect which is more evident in the high-mass 
tail. 

The other estimates, produced by Hyperzmass, are in quite 
good agreement with each other. The mass functions are slightly 
higher (on average by about 0.1 dex) when obtained using SED 
fitting procedures that can choose between two values of metal- 
licity. In fact, in this case, red galaxies can be fit with 0.2 Z and 
older ages, consequently resulting in larger stellar mass values. 
The effect of the extinction law is instead marginal. 



Article number, page 10 of 18 



I, Davidzon et al.: VIPERS: galaxy stellar mass functions 



-3 

,_, -4 

i 

X 
a) 

-5 



-5 
-6 



»o B „ 0.5<z< 0.7 
| 


. 0.7<z< 0.9 
w V 

i V 


O SB, Z sun V 
□ SB, Z, UI ,+0.2Z, un q 
SB + SMC, z mn <& 
A SB + SMC, Z sun + 0.2Z sun 
- V MAGPHYS 


A 

V 


i 1 i i i i 1 i i i i 1 i i i i 1 i i i i 1 

0.9<z< 1.1 


. ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 

l.Kz< 1.3 


Q 




o v 

V 

□ 

A 


A 

<^ 

a 

o 

v°v 



10.0 10.5 11.0 11.5 12.0 10.0 10.5 11.0 11.5 12.0 



log(Jl) [h 70 ~ 2 Jl e ] 

Fig. 7. Dependence of the mass function on the details of the stellar 
mass estimates, considering five different cases. Specifically, the points 
correspond to different choices of the stellar population synthesis code, 
metallicity (Z), extinction law (SB+SMC or SB alone) or the addition 
of secondary bursts to the smooth star-formation histories. Four cases 
correspond to SED fitting using Hyperzmass, for which the values of 
the adopted parameters are given at the bottom-left of the first panel. 
For details about the parameters adopted for MAGPHYS (downward tri- 
angle), we refer to Sect. 2.3. 



4. Comparison to previous work 

In this section we compare the VIPERS GSMF with other mass 
functions derived from different galaxy surveys (Sect. 4.1) and 
various semi-analytical models (Sect. 4.2). In the present paper 
we only exploit a sample that is approximately 60% of the final 
VIPERS catalogue, but this still provides results of a new level 
of accuracy, especially for the most massive galaxies. 

4.1. Comparison with observations 

We compare here our estimate of the GSMF with results from 
other galaxy surveys. We correct GSMFs (if necessary) to be in 
the same cosmological model with Q,„ = 0.3, O,,, = 0.7, h = 
0.70, and Chabrier (2003) IMF. We modify our redshift binning 
to be similar to other work. 

We choose from the literature eight surveys that adopt com- 
parable z-bins, half based on photometric redshifts (Fontana 
et al. 2006; Perez-Gonzalez et al. 2008; Ilbert et al. 2010; Bielby 
et al. 2012) and half on spectroscopic redshifts (Fontana et al. 
2004; Pozzetti et al. 2007, 2010; Moustakas et al. 2013). The 
spectroscopic redshift sample used by Moustakas et al. (2013) is 
obtained through a pioneering technique based on a low disper- 
sion prism and slitmasks (see Coil et al. 201 1), which results in 
a precision of cr. - 0.007(1 + z) (for their high quality sample 
(2^3, see Cool et al. 2013), i.e. comparable to the precision 
of the best photometric redshifts (Ilbert et al. 2013, who obtain 
cr : = 0.008(1 + z) and a very small percentage of outliers). 



The redshift ranges of the GSMFs shown in Fig. 8 are 0.4 < 
z < 0.6, 0.6 < z < 0.8, 0.8 < z < 1.0, with the exception of 
PRIMUS (Moustakas et al. 2013), which is at 0.5 < z < 0.65, 

0. 65 < z < 0.8, 0.8 < z < 1.0, and the first bin of VIPERS 
(i.e., 0.5 < z < 0.6). In the case of Bielby et al. (2012), who 
provide the GSMFs in four CFHTLS-Deep quadrants, we do not 
take into account the results in the Dl field, which is included 
in Wl field of VIPERS, nor from the D2 field, since in that re- 
gion we already have two independent estimates (i.e. COSMOS 
and zCOSMOS Ilbert et al. 2010; Pozzetti et al. 2010). For the 
VIPERS GSMFs we plot error bars accounting only for cr Po i ss , 

1. e. without adding the uncertainty due to cosmic variance, in or- 
der to be consistent with most of the literature data, for which 
only Poisson errors are available. 

Our results are in reasonably good agreement with most of 
the other GSMFs, with the exception of Moustakas et al. (2013), 
as visible in Fig. 8. 

At 0.8 < z < 1.0, a difference with Ilbert et al. (2010) and 
Pozzetti et al. (2010) is noteworthy: the likely reason is the pres- 
ence of a large structure, detected in the COSMOS field (Kovac 
et al. 2010), demonstrating the importance of the cosmic vari- 
ance in this kind of comparison. 

Some discrepancy remains with the estimates by Moustakas 
et al. (2013), whose mass functions appear to be systematically 
higher than all the previous measurements. A possible expla- 
nation could be related to the difficulty of estimating statistical 
weights for the PRIMUS survey: at magnitudes i =; 22.5 the SSR 
of PRIMUS is approximately 45% and at the limit of the survey 
it drops below 20% (Cool et al. 2013). Instead, in VIPERS the 
SSR is ~ 75% up to our magnitude limit i = 22.5 and to z — 1, 
making the statistical weight corrections small and very reliable. 
It is unlikely that the disagreement could be ascribed to cosmic 
variance, since PRIMUS surveyed a large area, with a number 
of independent fields (five fields with a total of 5.5 deg 2 ), which 
should significantly reduce this problem. 

Another source of difference between the observed GSMF 
in Fig. 8 is the choice of SEDs. Most work has used templates 
based on the BC03 library (Fontana et al. 2004, 2006; Pozzetti 
et al. 2007, 2010; Ilbert et al. 2010; Bielby et al. 2012), while 
Moustakas et al. (2013) derived their reference SEDs accord- 
ing to the stellar synthesis population model by Conroy & Gunn 
(2010). Perez-Gonzalez et al. (2008) also used a different library, 
selecting their templates from the PEGASE stellar population 
synthesis model (Fioc & Rocca-Volmerange 1997) and used a 
training set of ~ 2000 galaxies with spectroscopic z and a wide 
photometric baseline. 

In Fig. 8, the high statistical accuracy of VIPERS can be ap- 
preciated with respect to previous spectroscopic surveys. In par- 
ticular, VIPERS data provide tight constraints on the high-mass 
end of the GSMF. In fact, previous spectroscopic surveys like 
K20 and VVDS-Deep (Fontana et al. 2004; Pozzetti et al. 2007, 
respectively) were unable to probe this portion of the GSMF 
(log(A / t/A / l ) > 11.5) because of their relatively small area (52 
and 1750arcmin 2 respectively). Instead, GSMFs derived from 
photometric redshift surveys are characterised by a Poisson noise 
that is in general comparable to the level in VIPERS (Perez- 
Gonzalez et al. 2008; Ilbert et al. 2010), but they can be affected 
by the failures on photometric redshift estimates: even a small 
fraction of catastrophic redshift measurements can be relevant 
at large masses (Marchesini et al. 2009, 2010). Moreover, the 
sky area generally covered by high-z photometric surveys is not 
large enough for cosmic variance to be negligible. 

We will postpone a detailed analysis of the evolution of the 
local GSMF to a future work: differences in the details of the 
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Fig. 8. The VIPERS galaxy stellar mass functions from z = 0.5 to 1 (filled red circles, with a red shaded area accounting for the Poisson 
uncertainty). The 1/V max determinations of previous surveys are also shown by different symbols, along with their respective Poisson error bars. 
In the left-hand panel, whereas the VIPERS range is 0.5 < z < 0.6, the other GSMFs are computed between z = 0.4 and 0.6, with the exception of 
Moustakas et al. (2013) for which is 0.5 < z < 0.65, 0.65 < z < 0.8, 0.8 < z < 1.0; notice the very small error bars of the VIPERS data, despite 
the narrower redshift range. In the other two panels the bins of VIPERS are the same of the other surveys; also at these higher redshifts the error 
bars of the VIPERS GSMF are outstanding compared to them. 



available estimates from 2dFGRS, SDSS and GAMA (see Cole 
et al. 2001; Bell et al. 2003; Panteret al. 2004; Baldry et al. 2008; 
Li & White 2009; Baldry et al. 2012, providing quite differentre- 
sults especially at the massive end) are preventing a robust com- 
parison with our data. Only computing stellar masses and mass 
functions in a self-consistent way can provide constraints on the 
evolution of the GSMF down to z = (Moustakas et al. 2013). 

4.2. Testing models 

Besides the comparison with other surveys, it is important to 
check the agreement of our results with simulations. In this pa- 
per we limit ourselves to a preliminary analysis. Nevertheless, 
this first test provides intriguing results. A future work will be 
dedicated to more detailed investigations. 

The semi-analytical models (SAMs) we consider herein rely 
on the halo-merger trees of the Millennium Simulation (MS 
Springel et al. 2005) and the Millennium-II Simulation (MSII 
Boylan-Kolchin et al. 2009); namely, the MS (comoving box 
sizeL = 714 Mpc/i^J, particle mass = 1.23xl0 9 At Q hzl) is used 
by three of them (Bower et al. 2006; De Lucia & Blaizot 2007; 
Mutch et al. 2013) while the last one (Guo et al. 201 1) is based on 
the MSII (L = 143 Mpc , particle mass = 9.83 x 10 6 M Q h^). 
The tight constraints posed by VIPERS can be very useful when 
studying whether these models adequately reproduce the real 
universe. 

In Fig. 9, we show the mass functions derived from the mod- 
els of Bower et al. (2006); De Lucia & Blaizot (2007); Guo et al. 
(201 1), together with the VIPERS results. All the model GSMFs 
are computed from snapshots at the same redshifts. The narrow 
redshift binning we can do in VIPERS (Az = 0.1) allows us to 
compare simulated galaxies to observed ones at cosmic times 
very close to each other. For all of the three SAMs, we observe 
that the low-mass end of the GSMF is over-estimated. Low- 
mass galaxies seem to be in place too early. Such a discrepancy, 
already observed in other work (Somerville et al. 2008; Cira- 
suolo et al. 2010), is mainly due to an over-predicted fraction 
of satellite galaxies at those mass scales. This can be caused 
by an under-efficient supernova feedback and/or some issue in 
how the star formation efficiency is parametrised at high red- 



shifts (Fontanot et al. 2009; Guo et al. 201 1). Rescaling the sim- 
ulations to an up-to-date value of o"8 (in MS it is equal to 0.9), 
with the consequence of reducing the small-scale clustering of 
dark matter haloes, alleviates the tension only in part (Guo et al. 
2013). 

At higher masses, the VIPERS data provide very tight con- 
straints. At a first glance, De Lucia & Blaizot (2007) and Bower 
et al. (2006) seem to be in agreement with the observed GSMFs, 
while the Guo et al. (201 1) mass function lies systematically be- 
low by ^ 0.4 dex. However, it should be emphasised that those 
SAMs do not mimic the observational uncertainties on stellar 
mass. Adding this kind of error would increase the density of 
massive objects in the exponential tail of the mass function, and 
therefore the De Lucia & Blaizot (2007) and Bower et al. (2006) 
results should be considered at variance with observations also 
at log(A / t/A / l Q ) > 11. For this reason we do not show the De Lu- 
cia & Blaizot (2007) and Bower et al. (2006) GSMFs convolved 
with the observational uncertainty on the stellar mass. 

The effect of introducing observational uncertainties in the 
models is shown in Fig. 9 only for the Guo et al. (2011) mod- 
els, which foresee a lower density of objects in the massive end. 
We recomputed the Guo et al. GSMFs after convolving stel- 
lar masses with a Gaussian of dispersion 0.15 dex. The predic- 
tions of Guo et al. (2011) become then in fair agreement with 
VIPERS. With respect to De Lucia & Blaizot (2007), the main 
distinguishing features of Guo et al. (2011) are the high effi- 
ciency of supernova feedback and a lower rate of gas recycling 
at low mass. The transition from central to satellite status in the 
Guo et al. prescription also differs, resulting in a larger number 
of satellite galaxies than in De Lucia & Blaizot model. 

In several recent studies, the parameters of the models we 
discussed are tuned at z — by means of Bayesian approaches, 
fitting their stellar or luminosity function with data (Henriques 
et al. 2009; Bower et al. 2010; Lu et al. 2012). Observation at 
z > are then used to verify the model predictions (e.g. Lu 
et al. 2012). A different kind of calibration has been proposed 
by Mutch et al. (2013), who modify the input parameters in the 
SAM of Croton et al. (2006) to match observations at z = 
and z - 0.8 simultaneously. The VIPERS data will be crucial 
for such a test. Moreover, thanks to the narrow redshift bins of 
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Fig. 9. Comparison of the VIPERS mass function (red points, as in 
Fig. 4) with the semi-analytical models of De Lucia & Blaizot (2007); 
Bower et al. (2006); Guo et al. (2011) (grey solid, blue dashed, green 
short-dashed lines), whose GSMFs have been derived directly from the 
tables available in the Millennium database (Lemson & Virgo Consor- 
tium 2006). The Guo et al. (2011) stellar masses have also been con- 
volved with a Gaussian of dispersion 0.15 dex, to reproduce observa- 
tional uncertainty on stellar mass determinations; the resulting GSMFs 
are represented with green dot-dashed lines. 



VIPERS, one could precisely characterise the modelled GSMF 
between z = 0.5 and 1, by constraining the SAM parameters 
in small step of Az. That method could shed light on the time 
scale of the physical mechanisms that determine the evolution at 
higher masses (for instance, the AGN-feedback efficiency). 

The results obtained by Mutch et al. (2013) are compared to 
the VIPERS mass functions in Fig. 10. The plot shows a rea- 
sonable agreement beyond M - 10 11 A1 Q , not only at the red- 
shift of calibration (z - 0.83) but also in the other bins. The 
authors do not convolve their mass functions with a Gaussian 
uncertainty on stellar masses, because at least part of the uncer- 
tainties this procedure account for should already be included in 
the observational constraints they use. The Mutch et al. (2013) 
model is calibrated at z = 0.83 by using the results of Pozzetti 
et al. (2007), Drory et al. (2009), and Ilbert et al. (2010). Among 
these three GSMFs, only Pozzetti et al. (2007) is based on spec- 
troscopic data (VVDS-Deep), which are unfortunately quite lim- 
ited at higher masses. The other two estimates (Drory et al. 
2009; Ilbert et al. 2010) are derived from the COSMOS survey, 
which contains a significant over-density at z - 0.8. The strategy 
adopted by Mutch et al. to combine such information may lead to 
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Fig. 10. Comparison of the VIPERS mass function (red points) with 
the semi-analytical model of Mutch et al. (2013) (green shaded area at 
95% confidence limits). In several redshift bins Mutch et al. GSMF 
does not reach masses as high as VIPERS because the volume of the 
simulation (with a comoving box size L = 89.3 Mpc Jul) is smaller. In 
the right-hand middle panel (0.8 < z < 0.9) it is also shown, with a grey 
dashed line, the mass function that Mutch et al. obtain by combining ob- 
servational data from three different surveys (Pozzetti et al. 2007; Drory 
et al. 2009; Ilbert et al. 2010, grey triangles,diamonds, and squares rep- 
resent respectively). In addition, the yellow shaded regions represent 
the dispersion of the mass functions derived from the 57 SHMR mocks 
(see Sect. 3.3), in the same redshift bins of the VIPERS ones. 



overconfidence in the adopted constraints especially in the high- 
est mass range, where observations are most difficult. The ob- 
servational constraints the authors place on the high-mass end of 
the model (shown in Fig. 10) could be biased to some extent. To 
reconcile SAM and observations at log(Al/A1o) > 10.8, Mutch 
et al. have assumed a star-formation efficiency much higher than 
the one imposed by Croton et al. (2006), resulting in some ten- 
sions concerning the physical prescriptions of the model. In- 
triguingly, we note that the authors would significantly relieve 
these tensions if they were to add VIPERS data to their analysis. 

From a different perspective, the SHMR mocks we intro- 
duced in Sect. 3.3 are also calibrated at multiple z- We de- 
cide to test their reliability by deriving their GSMFs (Fig. 10). 
The agreement is remarkable: VIPERS data confirm the validity 
of the stellar-to-halo mass relation of Moster et al. (2013) that 
was use to construct these mocks. This relation connects galax- 
ies with their hosting dark matter halo by means of a redshift - 
dependent parametrisation that has been calibrated through the 
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GSMFs of Perez-Gonzalez et al. (2008) and Santini et al. (2012) 
up to z — 4. Because of the lack of tight constraints used by 
Moster et al. for the most massive galaxies (the data from Perez- 
Gonzalez et al. 2008 have lower statistics than ours), the SHMR 
mass functions diverge at high mass from our estimates. 

5. Evolution of the mass function of the red and 
blue galaxy populations 

In order to distinguish the contribution of quiescent and actively 
star forming galaxies to the global evolution, we now split the 
sample according to the galaxy rest-frame (U - V) colour (see 
Fritz et al. (2013) for extensive discussion). 

5.1. Classification of galaxy types 

The absolute magnitudes for galaxies in the VIPERS catalogue 
were computed from the same SED fitting procedure described 
in Sect. 2.3, applying a k- and colour-correction, derived from 
the best fit SED, to the apparent magnitudes in the bands that 
more closely match the rest-frame emission in the U and V filters 
(see Fritz et al. 2013, for details). In this way, (U — V) rest-frame 
colours can be reliably computed within the redshift range of the 
survey, evidencing the classical bimodality (Fritz et al. 2013) and 
allowing us to separate red-sequence from blue-cloud galaxies 
(Strateva et al. 2001; Hogg et al. 2002; Bell et al. 2004). 

The valley between the two populations is found to be 
slightly evolving toward bluer colours at earlier epochs. This 
photometric classification, despite its simplicity, can be consid- 
ered as a good proxy for the selection of quiescent and star- 
forming galaxies, as discussed by Mignoli et al. (2009) using 
zCOSMOS data: 86% (93%) of the galaxies selected as being 
photometrically red (blue) are also quiescent (star-forming) as 
defined from their spectra. 

To verify and validate our selection method, we also derive 
galaxy photometric types by fitting our photometry with the em- 
pirical set of 62 templates used in Ilbert et al. (2006), which 
was optimised as to to refine the match between photometric and 
spectroscopic redshifts in the VVDS. The same set was also used 
to classify galaxies in several otherpapers (e.g. Zucca et al. 2006, 
2009; Pozzetti et al. 2010; Moresco et al. 2010). The classi- 
fication of VIPERS galaxies resulting from this second method, 
matches reasonably well with the (U-V) colour selection. More 
than 70% of the red galaxies are defined as early-type objects by 
the SED analysis, while more than 95% of blue galaxies are clas- 
sified as late-types. For red galaxies this worsen beyond z = 1.1, 
where only 55% of the red galaxies are classified as early-type 
in terms of their SED. In the same redshift range, instead, 98% 
of blue galaxies are classified as late-type objects. 

5.2. Blue and red galaxy stellar mass functions 

Using this classification, we are now in the position to quantify 
the contribution of red and blue galaxies to the GSMF and in 
particular to its high-mass end. The results are shown in Fig. 1 1 . 
The mass functions for each class are estimated in bins of 0.2 in 
log(Al Q ), using the same 1/V ma x method described in Sect. 3.2. 
Fits with the usual Schechter function are provided, as described 
in the caption, to highlight evolution (or absence thereof) as a 
function of redshift. 

The predominance of red (passive) objects among the mas- 
sive galaxies is clearly visible in all redshift bins, with blue (ac- 
tive) galaxies mainly contributing at lower masses (Ai < Al*). 



Since the mass completeness limit Alum for the blue popula- 
tion extends to sufficiently small masses, we can perform the 
Schechter fit letting Al*, <£>*, and a free. The slope of the low- 
mass end remains almost constant in redshift for the blue pop- 
ulation, with 1.2 < a < 1.3, up to z - 0.9, as seen in previous 
works (cf. Pozzetti et al. 2010). At redshift larger than this it 
can no longer be constrained. Because of the high values of 
the mass completeness limit AW for the red population (see 
Sect. 3.1), it is not possible to study the red sample in the same 
range, e.g. determining the evolution of a (Ilbert et al. 2010) or 
a possible upturn of the GSMF (cf. Drory et al. 2009). 

From these measurements we can determine the value of 
Across, where the blue and red GSMFs intersect, i.e. the divid- 
ing line between the ranges in which blue and red galaxies re- 
spectively dominate the mass function (Kauffmann et al. 2003b). 
The physical meaning of Across has been questioned (Bell et al. 
2007), but it is in general considered as a proxy to the transi- 
tion mass of physical processes such the quenching of star for- 
mation, (responsible for the migration from the blue cloud to 
the red sequence), or the AGN activity (e.g. Kauffmann et al. 
2003a). Moreover, its clear dependence on environment (Bol- 
zonella et al. 2010) points to an interpretation of the galaxy trans- 
formation linked not only to secular processes. 

We quantify the value of the transition mass in each red- 
shift bin using the 1 / V ma x measurements. The transition mass 
increases from log(Alcross/At Q ) = 10.4 at z ^ 0.55 to 
log(At C ross/Alo) = 10.6 at z 0.75, as shown in Fig. 12. This 
trend is very well fitted by a power law oc (1 + z) 3 . Beyond 
z — 0.8 our Across estimates (open circles) should be formally 
considered as upper limits, since they fall below the mass com- 
pleteness limit of red galaxies, but at least up to z = 1.0 they can 
be considered as a good approximation of the real values, given 
their proximity to the limit. 

In Fig. 12 we also plot results from previous studies. In this 
respect, it is important to underline that the value of the com- 
pleteness mass Alcross provided by the various authors can differ 
significantly from each other, depending on the adopted classi- 
fication. For instance, the results of the morphological classifi- 
cation used by Bundy et al. on the DEEP2 survey fall outside 
of the mass ranges considered in the plot. This could be related 
to part of the "red and dead" galaxies at such redshifts becom- 
ing ellipticals (in morphological sense) at a later stage (Bundy 
et al. 2010). In fact, when we split the DEEP2 sample on the 
basis of the (U - B) bimodality (upper half-circles), the results 
are then in agreement with our findings. Our estimates of Alcross 
are fairly consistent with those of Vergani et al. (2008), Pozzetti 
et al. (2010), Moustakas et al. (2013). The estimates by Vergani 
et al. (2008) (astriskss in Fig. 12), rely on the identification of 
the D4000 break in the VVDS spectra and have a steeper red- 
shift evolution, At cro ss K (1 + z) 4 - Pozzetti et al. (2010) derived 
A( cr oss from the GSMFs of the zCOSMOS (lOk-bright) sample 
split using different criteria: a cut in specific SFR (i.e. sSFR 
= SFR/At % 10 _1 Gyr _1 - open squares in the plot), morphology 
(spheroidal vs disc/irregular galaxies - diamonds), and best-fit 
SEDs (same photometric types discussed in Sect. 5.1 - trian- 
gles). Moustakas et al. (2013) define star-forming galaxies as 
lying in the so-called main sequence of the SFR (estimated from 
the SED fitting) vs Al diagram (Noeske et al. 2007). They find a 
flatter evolution, with Al cro ss K (1 +z) 15 (dashed line). 

5.3. Evolution of the blue and red populations 

To study the processes of star formation quenching that cause 
the transition of galaxies from the so-called blue cloud to the 
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Fig. 11. The galaxy stellar mass functions of the blue and red populations in VIPERS, derived using the 1/K max . Symbols (circles and diamonds 
respectively) are filled for data above the corresponding completeness limit M\, m (vertical lines) and empty below. Error bars account for Poisson 
noise only. The Schechter fit of the two populations in the bin 0.5 < z < 0.6 (solid blue and red lines) is reported for reference as a dashed line in 
the other panels. The solid black line in each panel give the Schechter best-fit to the whole VIPERS sample in that redshift bin. 



red sequence (Faber et al. 2007), we measured the evolution of 
the galaxy number density of blue and red populations, namely 
p^J ue (z) and p^ d (z). These estimates are derived using the 1 / V max 
method, taking into account both Poisson noise and cosmic vari- 
ance. We also verified, however, that the results would be essen- 
tially the same if we had measured number densities by integrat- 
ing the Schechter best-fitting functions. We explore four narrow 
bins of stellar mass to highlight the dependence of the quench- 
ing processes on this parameter. To improve statistics at large 
stellar masses, here we choose wider redshift bins: 0.5-0.7, 0.7- 
0.9, 0.9-1.1, 1.1-1.3 (the same used by Fritz et al. 2013, for the 
analysis of the VIPERS luminosity function). 

At intermediate masses (10.8 < log(A1/A1o) < 11.1) the 
number density of red galaxies p^ d increases by a factor of 
~ 2.5 from z - 1 to z — 0.6, whereas at higher masses 
the variation is much smaller: red galaxies with mass 11.1 < 
log(A1/ M Q ) < 11.4 evolve from a number density p^ d = 

(5.3 ± 0.6) x 10~ 5 Mpc~ 3 h 3 10 to (9.5 + 1.1) x 10~ 5 Mpc~ 3 h] Q in 
the same redshift interval (about 80% increase). This result is in 
line with the mass-assembly downsizing scenario highlighted in 
previous works (Cimatti et al. 2006; Pozzetti et al. 2010; Ilbert 
et al. 2010): elliptical galaxies with M > 10 11 M@ built their 
stellar mass well before the less massive ones. With the VIPERS 
data we are able for the first time to provide significant evidence 
of this trend for such massive galaxies log(A1/A1o) > 11.4 at 
these redshifts. Thanks to the combination of large volume and 
sampling of VIPERS we can look for evolutionary trends in this 
mass range, whithout being affected by the cosmic variance. We 
find that the number density of the most massive red galaxies 



in our redshift range increases by 45%, confirming that they did 
not experience a strong evolution between z — 1.2 and z — 0.6. 
Quenching mechanisms seem to be more efficient at low and 
intermediate masses, as recently suggested also by Moustakas 
et al. (2013). With respect to PRIMUS, the VIPERS survey ex- 
tends this finding to higher masses (^(AI/AIq) > 11.4) and 
redshifts (up to a* 1.2), The evidence of mass dependence of 
quenching is in agreement with Peng et al. (2010), although 
other mechanisms could play a non-negligible role (e.g. galaxy 
mergers, Xu et al. 2012). 

The co-moving number density of blue galaxies, instead, is 
found to be relatively stable between z - 1 to z — 0.6 for ob- 
jects with mass 10.5 ^ log(A4/A / l ) < 10.8, with a 10% vari- 
ation. For larger-mass blue galaxies, for which the sample is 
complete at all redshifts, the density p^J ue of objects with mass 
10.8 < log(A1/A1o) < 11.1 shows indications for a mild in- 
crease between z — 1.2 and z = 0.8. 

The most massive blue galaxies (log(A / t/A / l ) > 11.4) dis- 
appear downward of z - 0.6 (see the right panel in Fig. 13), 
giving us a hins about the time at which the stellar mass cut-off 
of active galaxies developed. When the whole VIPERS sample 
will be available, we will carry on the analysis of the massive 
end buildup with more robust statistics. Moreover, a step for- 
ward for a better comprehension of this picture will be the use of 
spectral features to determine reliable estimates of the SFR, and 
therefore to better separate passive and active galaxies. 
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Fig. 13. Evolution of the abundance (number density) of the blue and red galaxy populations with different stellar masses in VIPERS (filled 
circles and diamonds respectively). The error corridors reflect the overall uncertainties, which include both Poisson noise and cosmic variance 
added in quadrature. 
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Fig. 12. The values of the transition mass Across as computed from 
Fig. 11, plotted as a function of redshift. The VIPERS measurements 
are given as black open circles, with a downward arrow depending on 
whether the transition mass is above or below the completeness mass 
of at least one of the two classes. The solid line is a fit with a (1 + z) 3 
power-law to the VIPERS points between z = 0.5 and z = 0.8. These 
are compared to literature estimates in grey. Points from Pozzetti et al. 
(2010) are obtained using three different classifications: a separation ac- 
cording to specific SFR (diamonds), a best-fit SED classification (trian- 
gles), and a morphological classification (squares). The points of Bundy 
et al. (2006) are based on the (U-B) bimodality and [Oil] emission re- 
spectively (see text). The points by Vergani et al. (2008) are based on a 
spectral classification (D4000 break). The value from PRIMUS (Mous- 
takas et al. 2013) at z = 0.9 is reported as a cross, while the dashed line 
traces an evolution oc (1 + z) , as suggested in that paper; these authors 
classified active and quiescent galaxies with respect to their position in 
the SFR vs M diagram. 



6. Conclusions 

We have measured the GSMF between z = 0.5 and z = 1 .3 us- 
ing the first data release of VIPERS. The forthcoming VIPERS 



Public Data Release 1 (PDR-1) will contain the catalogue of the 
53 608 spectroscopic galaxy redshifts used in the present anal- 
ysis. The galaxy stellar masses have been estimated through 
an SED fitting technique, relying on a large photometric base- 
line, and in particular on a nearly full coverage of our fields 
with near-infrared data; we performed several tests to verify that 
the systematics intrinsic to the method of SED fitting (e.g. the 
parametrisation of the SFH) are not introducing any bias in our 
analysis. The large volume probed by VIPERS results in ex- 
tremely high statistics, dramatically reducing the uncertainties 
due to Poisson noise (cr Po ; ss ) and sample variance (cr cv ). We es- 
timate the latter by using 57 galaxy mock catalogues based on 
the MultiDark simulation (Prada et al. 2012) and the stellar-to- 
halo mass relation of Moster et al. (2013). These mocks closely 
reproduce the characteristics of the VIPERS survey. 

We empirically determined a completeness threshold Ainm 
above which the mass function can be considered complete: 
this limiting mass evolves as a function of z, ranging from 
log(A1/A1 ) = 9.8 to 11 in the redshift interval 0.5-1.1. We 
focused our analysis on the high-mass end of the GSMF, where 
VIPERS detects a particularly high number of rare massive 
galaxies. The main results we obtain are: 

- VIPERS data tightly constrain the exponential tail of the 
Schechter function, which does not show significant evolu- 
tion at high masses below z = 1.1. The same result is pro- 
vided by the analysis of the co-moving number density p^, 
calculated in different bins of stellar mass: at z — 1 .2 most of 
the massive galaxies with log(A1/A1 ) > 1 1.4 are already in 
place, whereas below k^Wl/Alg,) = 11.4 the galaxy num- 
ber density increases by a factor of ~ 3.5 from z — 1.2 to 
z * 0.6. 

- We compare our observed GSMFs with those derived from 
semi-analytical models (De Lucia & Blaizot 2007; Bower 
et al. 2006; Guo et al. 2011). While the discrepancy at low 
masses between models and observations is well established 
and has been exhaustively discussed in literature, predictions 
at the high-mass end of the GSMF have not yet been verified 
with sufficient precision. We show that the high accuracy of 
the VIPERS mass functions make them suitable for this kind 
of tests. From a first analysis, the VIPERS data appear to be 
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consistent with the Guo et al. (201 1) model, once the uncer- 
tainties in the stellar mass estimates are taken into account. 
A more detailed analysis will be the subject of a future work. 
We suggest that VIPERS GSMFs can be effectively used to 
constrain models at multiple redshifts simultaneously, as pi- 
oneered in the work of Mutch et al. (2013). 

- We divide the VIPERS sample by means of a colour criterion 
based on the (U - V) bimodality (Fritz et al. 2013) and we 
estimate the blue and red GSMF in the same range, 0.5 < z < 
1 .3. We find that the transition mass above which the GSMF 
is dominated by red galaxies is about 10 10 4 Ai at z — 0.55 
and evolves proportional to (1 + z) 3 . 

- The number density of the red sample shows an evolution 
that depends on stellar mass, being steeper at lower masses. 
At large stellar masses the quenching of active galaxies has 
not been thoroughly studied, because of the poor sample 
statistics. We obtain a first impressive result with VIPERS, 
detecting at z - la significant number of very massive active 
galaxies with log(A / l/A / l ) > 11.4, which have all migrated 
onto the red sequence by z = 0.6, i.e. in about 2 Gyrs. 

The first data release of VIPERS has thus allowed us to study 
the evolution of the galaxy stellar mass function over an unprece- 
dented volume at redshifts z = 0.5 — 1.3. We emphasise the con- 
straining power of this dataset, particularly for the abundance of 
the most massive galaxies, both quiescent and star-forming. In 
forthcoming studies we will make full use of the growing sam- 
ple and of the measurement of spectral features, in order to in- 
vestigate the cosmic star formation history and confront galaxy 
formation models at high redshift. 
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